Occupation numbers from functional integral 
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Occupation numbers for non-relativistic interacting particles are discussed within a functional 
integral formulation. We concentrate on zero temperature, where the Bogoliubov theory breaks 
down for strong couplings as well as for low dimensional models. We find that the leading behavior 
of the occupation numbers for small momentum is governed by a quadratic time derivative in the 
inverse propagator that is not contained in the Bogoliubov theory. We propose to use a functional 
renormalization group equation for the occupation numbers in order to implement systematic non- 
perturbative extensions beyond the Bogoliubov theory. 



I. INTRODUCTION 

Occupation numbers are the central tool for the sta- 
tistical description of non-interacting particles or quasi- 
particles. Extending this concept to interacting particles 
allows for the description of many interesting phenomena, 
as condensates, the formation of gaps and pseudogaps or 
nontrivial scaling for the momentum dependence of the oc- 
cupation numbers. Most of these phenomena occur at low 
temperatures T and we will concentrate on this case. 

For non-relativistic bosons and small couplings the use 
of perturbation theory, within a Hamiltonian framework 
with operator valued fields, has led to the successful Bo- 
goliubov theory [![. This theory breaks down, however, 
in low dimensions and for strong couplings. In one di- 
mension (d = 1) and for T = 0, the occupation number 
n(q) for a mode with momentum q diverges proportional 
to the inverse momentum, n(q) ~ = \q\- The parti- 

cle density is the momentum integral over the occupation 
numbers and diverges logarithmically in this case. Further- 
more, the Bogoliubov theory is built on a nonzero conden- 
sate density or order parameter (f>Q. Such an order param- 
eter leads to spontaneous symmetry breaking of the global 
[/(l)-symmetry associated to the conserved particle num- 
ber. It is known that for d = 1 and T = no spontaneous 
symmetry breaking of a continuous symmetry is possible - 
the order parameter 4>q has to vanish in the infinite volume 
limit. Also for d = 2 the situation is problematic. Now one 
finds a nonzero </>o for T = and the Bogoliubov approxi- 
mation to the particle density shows no infrared divergence 
anymore. However, for any nonzero temperature the order 
parameter must again vanish due to the Mermin- Wagner 
theorem 0, M T > 0) = 0. This reveals again a prob- 
lem for the perturbative treatment where n(q) ~ 4>o/q for 
small q. For q > the occupation numbers should be con- 
tinuous functions of T and this indicates a contradiction 
between the Bogoliubov approximation and the Mermin- 
Wagner theorem. 

Going beyond leading order perturbation theory (Bogoli- 
ubov approximation) becomes quite intricate in the opera- 
tor formalism. A good alternative is the functional integral 
formulation which has already given a rather detailed pic- 
ture of Bose-Einstein condensation for interacting bosons, 
largely by use of the renormalization group [|[. In this 



paper we address the computation of occupation numbers 
within a functional integral approach. 

Already at the very basic level of the definition of n(q) 
one needs to be careful to subtract "counterterms" that di- 
verge in the ultraviolet limit. We work with a momentum 
cutoff q 2 < A 2 and show that a too naive definition of n(q) 
would lead to a divergent particle density for A — > oo. One 
needs to introduce counterterms that arise from the tran- 
sition from the Hamiltonian formalism to the functional 
integral. With this correct starting point the functional 
integral formalism is well suited to solve the infrared prob- 
lems of the Bogoliubov theory, as well as potential ultravi- 
olet problems of extensions of it. 

For T — we show how the infrared problems of the low 
dimensional Bogoliubov theory can be resolved by an in- 
vestigation of the frequency and momentum dependence of 
the full propagator. In particular, it is crucial that the fre- 
quency dependence of the inverse propagator is quadratic 
in the range of small q and not linear as for the Bo- 
goliubov theory. The use of the propagator computed in 
[J] removes the infrared problems. For q — ► the particle 
number behaves now ~ V~ x l 2 q~ 1JrT i instead of ~ ^o9~ 1 m 
the Bogoliubov approximation. Here V is the coefficient of 



the term 



in the inverse propagator, with u> the fre- 



quency. The anomalous dimension rj vanishes for d — 2,3, 
while it is in the range < rj < 1/2 for d = 1. This renders 
the particle number finite for all dimensions. Furthermore, 
the leading behavior of n(q) for q — ► becomes indepen- 
dent of 0o and is therefore not affected by the vanishing 
of 0o at T = (for d = 1) or T > (for d = 2). In 
three dimensions and for a small coupling the behav- 
ior n ~ V~ 1 / 2 q~ 1 is relevant only for very small q, whereas 
the validity of the Bogoliubov approximation n ~ <f>o j q 
extends to a large range of q. 

Extensions beyond the Bogoliubov theory have also to 
cope with ultraviolet problems. They occur if the fre- 
quency dependence of the propagator is not taken into ac- 
count properly. We propose here a renormalization group 
treatment based on the average action 

HQ- An in- 
frared cutoff k for the fluctuations is introduced in order 
to extrapolate between the microscopic action for k = A 
and the full quantum effective action for k — 0. Corre- 
spondingly, we study flow equations for /c-dependent oc- 
cupation numbers nk{q)- They yield the physical occupa- 
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tion numbers n(q) as the IR-cutoff is removed for k = 0. 
The renormalization group flow permits an improvement in 
the computation of n(q) since for every momentum range 
only the effective renormalized couplings in this momen- 
tum range matter. In particular, we use the renormaliza- 
tion group flow in order to show two properties. The first 
is the non-renormalization of the condensate density for 
non-relativistic bosons. It is given by the unrenormalized 
order parameter, n c = 0q. Second, this method is free 
of ultraviolet problems. Our methods can be extended to 
coupled systems of bosons and fermions (cf. sect. IXIII and 
appendix B) which are relevant for the BEC-BCS crossover 
for ultracold fermions. There the issue of the condensate 
density is more complicated @. A renormalization group 
flow for the total density in such systems has already been 
investigated in [9(. 



This paper is organized as follows. In sect. |TT] we dis- 
cuss the occupation numbers in a functional integral con- 
text. They are related to a go-integral over the propagator, 
with qo the Euclidean frequency (Matsubara frequency for 
T > 0). This section also introduces the (ultraviolet diver- 
gent) countertcrms. Sect. [HI] specializes to zero tempera- 
ture and sect. IIVI deals with the case where the frequency 
dependence of the inverse propagator is linear. This yields 
the Bogoliubov approximation and first extensions of it. In 
sect. [V]we extend the general form of the inverse propa- 
gator by taking terms quadratic in the frequency (or also 
higher orders) into account. We show that this solves the 
infrared problems of the Bogoliubov theory in low dimen- 
sions. We present differential equations for the momentum 
dependence of the occupation numbers and the particle 
density. In sect. I VII we turn to a heuristic renormalization 
group improvement. This helps for a proper understanding 
of the occupation numbers for large momentum. 



The following sections introduce an exact renormaliza- 
tion group equation for the occupation numbers. It moti- 
vates the simpler heuristic renormalization group improve- 
ments. We start in sect. IVIII by introducing a generating 
functional for the occupation numbers. This can be made 
scale dependent, and an exact flow equation is derived in 
sect. IVIIII We perform a non-perturbative approximation 
by quadratic truncation in sect. HXl The flow shows a mu- 
tual influence of occupation numbers for different momenta 
- it cannot be separately described for individual momen- 
tum modes. The corresponding density transfer function 
is discussed in sect. [X] In sect. IXII we apply the functional 
renormalization group equations to non-relativistic bosons 
for T = 0. Sect. IXIII discusses gaps and pseudogaps in the 
light of renormalization group improvement and we give a 
short summary of our main results in sect. IXIIII 



II. OCCUPATION NUMBERS AND 
RENORMALIZATION 

We start with the partition function for a nonrelativistic 
bosonic particle 



Z = J V X exp(-S\x\), 



Six] = J (x*8rX+2]jfVx*Vx-crx*x) + Sini. W 
The complex field x may be expressed by its Fourier modes 



X (x) = X (t, x) = J e^ X (r, q)=J e^ X (q), (2) 

with 



q = (qo,q), \ = {^y d J d d q , 

= (27r)( d+1 ) / dq a I d d q. (3) 



For nonzero temperature T the Euclidean time r parame- 
terizes a circle with circumference fl T = T _1 and the Mat- 
subara frequencies qo — 2^nT , n € 2, are discrete, with 
J = T^2. The second thermodynamic parameter is the 

n 

effective chemical potential a. We regularize the theory 
by a momentum cutoff q 2 < A 2 and take A — * oo when 
appropriate. 

Furthermore, we assume the invariance of the classical 
action S under a global abelian symmetry of phase rota- 
tions x e X) corresponding to a conserved total particle 
number 



N 



n(x) = fl d n(q) 



(4) 



with fid the volume of d-dimensional space (fld+i = Qd^r)- 
In our normalization n(q) is dimensionless. The momen- 
tum distribution n(q) denotes the fraction of particles with 
momentum q. This is the quantity we want to investigate 
in this paper. For definiteness, we may consider a torus 
with finite fid and discrete momenta, taking fid — > oo at 
the end. For discrete momenta we write 



(5) 



and identify n(q) with the occupation number in momen- 
tum space. 

Following the Noether construction we can express n(x) 

as 

= YT T Jdr(x*(r,x)x(r,x)) - \ J ±jL (6 ) 



and, correspondingly 
1 



n(q) = ^-j- 1 (x*(qo,q)x(qo,q)) - 2 - 



(7) 



3 



Note that n(q) is not the Fourier transform of n(x). 

The constant part —1/2 in n(q) arises from the connec- 
tion between the functional integral and the operator for- 
malism Q. For a single degree of freedom it is the sym- 
metric combination of annihilation and creation operators 
(a ji a + aa i )/2 that translates into x*X- Therefore the num- 
ber operator n — a) a = (a) a + aa^)/2 — 1/2 acquires an 
additional contribution and this generalizes to all momen- 
tum modes. 

In the following we will work with a basis of real fields 
Xi,X 2 defined by x(%) = ^%(xi(x) + iXi{x)), such that 
Xa(-q) = X*a(i) and 



n(q) 



1 1 



2 2fi(2+i 



(<x!(9)Xifa)> + <x2(9)Xa(«)>)- (8) 



The connected part of the two point function describes the 
propagator Q 

(x* a (<l)Xb(q r )) = Qab(q,q) + (x* a (q))(xb(q))- (9) 

For a translation invariant setting Q is diagonal in momen- 
tum space 



Gab(q, q') = G ab (q)S(q - q'), 



(10) 



with 5(q - q') = {2Tr) d+1 5{q - q' )d d (q- q'). Also, trans- 
lation invariance implies for a possible order parameter 
(Xa(q)) = V^4>o5{q)8 a i with real 0o- Here we have cho- 
sen the expectation value in the one-direction without loss 
of generality. We define 



m = ^(Gn(q)+G 22 (q)) 



(11) 



and express n(q) in terms of the full propagator G(q) and 
the order parameter <j> 



n(q) 



g(qo,q) +^>tKq) 



n P (q) + 4>l s (q)- 



(12) 



In the absence of spontaneous symmetry breaking (0o = 
0) there is no difference between xi and \i- Since the 
expectation values of the type (xx) an d (x*X*) vanish in 
this case, it is often convenient to work in the "complex 
basis" of a one component complex field. The propagator 
is then associated with 



(x(q')x(q))=G(q)S(q-q / ) 



(13) 



and G(g) will typically be a complex function. For free 
bosons one has 



q 



G (q) = \ i i° + ^jj- (T 



I G (q) = T^2(2TTinT + q 2 /2M-a)- 1 . (14) 



Performing explicitely the Matsubara sum yields 



i(q) = / G Q (q) - 



1 



exp 



( q 



\ 2MT T 



- 1 



(15) 



and we recover the usual occupation number for free 
bosons. Notice that the subtraction of the term —1/2 in 
eq. (fT2|) is crucial for the vanishing of n(q) for q 2 — > oo. 

For fermions the situation is similar. There is no 
fermionic expectation value and the additive constant 
changes sign according to the anticommuting properties 
^(a F aF — a-FOp) = np — i, i.e. 



n F {q) 



G F {qo,q) + 



(16) 



The Matsubara sums for T^O involve now half-integer n 
and we recover for a free theory the Fermi distribution 



n F (q)= exp 



q 



2MT 



1 



(17) 



In presence of interactions, however, the exact propa- 
gator G{q) can be a rather complicated object. We dis- 
tinguish between the vacuum and "occupied states" where 
n(q) differs from zero. Often, but not always, an occupied 
state will go along with a condensate 0- F° r non- 

relativistic interacting particles the propagator becomes 
complicated only for the occupied states. The complexity 
of the occupied states is due to the presence of a new scale 
related to the density. In contrast, the vacuum of nonrel- 
ativistic theories is characterized by a simple propagator 
fill , 0] , even in the presence of interactions. For the vac- 
uum, we will use a general non-renormalization property of 
the wave function renormalization in order to clarify some 
relations between the microscopic formulation of the theory 
and macroscopic properties. 

Let us consider T = where the qo integration is continu- 
ous. (More precisely, we define the result of a go-integral as 
the principal value and J qo f(q ) = ^im^j ^ f^/(?o)-) 

The vacuum is defined by the limit n{q) — > 0. Since no 
condensate exists in vacuum (<f>g — 0) we will work in the 
complex basis. For a stable particle and q 2 > the full 
vacuum propagator G(q) must have a pole for purely imag- 
inary qo with positive definite imaginary part. (For free 
particles the location of the pole is qo = iq 2 /(2M) for 
a = 0.) Often we can write the go _m tegration as 



= 5r / d qo = T-0 

qo Z7r J ^ ./upper 

— oo 



) 



(18) 



lower 



where f u denotes the integral closed on the upper half 
plane in complex go-space. As long as the pole remains in 
the upper half plane the result of the go integration is given 
by the residuum of the pole 



/ G{q^q)= l -ZZ 1 {q). 
Jqo z 



(19) 



Here we define the wave function renormalization Z^, (q) by 
the derivative of the inverse propagator at the location of 
its zero 



d 



G _1 (w (g),g) = 0, to = -iq . 



(20) 
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For a pole in the lower half plane the r.h.s. of eq. (fl9|) 
picks up a minus sign. 

Eq. (|19p constitutes an important relation between the 
occupation number and the behavior of the full propagator 
at its pole 

n{q) = \(Z-\q)-l). (21) 

From eq. (|2ip we can derive immediately a crucial non- 
renormalization property for the vacuum of a theory for 
nonrelativistic particles. Since for the vacuum n(q) must 
vanish one concludes that Z w (q) is not renormalized. It 
keeps its classical value even in presence of interactions 

lim Z u (q) = 1. (22) 

This non-renormalization property resolves several puz- 
zles that would arise if the interaction effects would lead to 
a nontrivial renormalization ^ 1. The first concerns the 
normalization of the particle number. Indeed, the global 
[/(l)-symmetry implies a conserved total particle number 
N, but it does not fix the normalization of N. Obviously, 
ZN is also conserved for arbitrary constant Z. The ques- 
tion of the proper normalization of N can be addressed 
from a microscopic and a macroscopic point of view. In the 
operator language the normalization of the particle num- 
ber follows directly form the normalization of the creation 
and annihilation operators via the commutation relation 
[a, a^] — 1. Constructing the functional integral results in 
a normalization of the particle number through the normal- 
ization of the term x*d T X m the action. (For an arbitrary 
normalization of the fields the term linear in the frequency 
will take the form iZ\(q)qox* (q)x(l) instead of the nor- 
malization Za(<i) = 1 employed in eq. |T]). The particle 
number reflects then this normalization and the r.h.s of eq. 
([7]) is multiplied by a factor Z^{q).) This microscopic point 
of view relates the particle number to the full propagator 
G(q) for the bare or microscopic fields. 

On the other hand, from a macroscopic point of view it 
should be possible to express all physical quantities relat- 
ing to macroscopic scales in terms of the n-point functions 
for renormalized fields, without explicit reference to the 
microscopic formulation. This is a basic concept in quan- 
tum field theory and leads to the concept of renormalizable 
theories when the characteristic macroscopic length scale 
is many orders of magnitude larger than the microscopic 
length A -1 . In relativistic quantum field theories the renor- 
malized fields xr — Z X ^ 2 X often involve a wave function 
renormalization Z than diverges for A — > oo. Our defini- 
tion of the particle number seems at first sight to be at 
odd with this general property of quantum field theory. It 
is expressed in terms of the bare propagator G(q) and not 
in terms of the two point function for the renormalized field 
which is ZG(q). Therefore the particle number seems to 
keep some memory of the microscopic normalization. 

In view of eq. (|22|) the resolution of the problem 
is simple: in vacuum the renormalized fields Xr(q) = 
^a>(?) x(<?) are identical to the microscopic fields. Only 
for nonzero T and n one will find Zu{q) ^ 1. However, this 



quantity can now be expressed in terms of macroscopic 
quantities which describe the difference between a thermo- 
dynamic equilibrium state and the vacuum. As it should 
be, no microscopic physics is involved. 

A second puzzle is related to a possible dependence of 
the microscopic wave function renormalization Z\(q) on 
the scale A (in some given arbitrary normalization of the 
field x)- Indeed, the definition of the functional integral 
does not only involve the action S but also a regularization 
scale A. Changing A results in a change of S such that for 
a given model the action S\ depends on the choice of A. 
If for a given normalization of the fields the wave function 
Za(q) would run with A, it would not be clear which Z\(q) 
should be used for the definition of n. Fortunately, the 
non-renormalization property in vacuum (|22[) implies that 
Zuj,a(q) = 1 does not depend on A. 

We close this section by a formal side remark concerning 
the qo-hrtegral of the vacuum propagator. For large real qo 
the renormalization effects vanish such that 

lim G- 1 (q ,q)=iq . (23) 

|go|^oo 

Using the identity (for e > 0) 

\= [ (m + c)- 1 (24) 

L Jqo 

we can write the go-integral in a manifestly convergent form 

[ G{q Q ,q)-\= ( {Giqo^-itqo + e)- 1 }. (25) 
Jqo Jqo 



III. OCCUPATION NUMBERS AT ZERO 
TEMPERATURE 

For T — and n(q) > we may employ eq. (|2"Tj) (or 
generalizations thereof) in order to relate the occupation 
number to the wave function renormalization Z^(q). Be- 
fore exploiting this relation in sects. IIXI and 1X1 for bosons, 
we first discuss the analogous case of fermions. In the ab- 
sence of spontaneous symmetry breaking one has a single 
complex propagator Gp- For a pole in the upper half plane 
the go-integration yields 

^=l(i- z ; 1 m> z ^= I z^$' (26) 

whereas for a pole in the lower half plane one finds 

n ^-V 1 + Z ^ Z ^-^^l (27) 

For a free theory and a > the pole remains in the upper 
half plane for q 2 >k F , k F = 2Ma. For these momenta 
n(q) = implies Z u {q) = 1. In contrast, for q 2 < k F the 
pole is in the lower half plane. Now n(q) = 1 results again 
in Z u {q) = 1. Thus a sharp Fermi surface is equivalent to 
Zu(q) = 1- In contrast, the presence of interactions can 
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"smoothen" the Fermi surface and lead to Z u (q) ^ 1, even 
for T = 0. 

The jump in the particle number np(q) is an interest- 
ing example how quantum mechanical discreteness arises 
from a perfectly continuous functional integral. Let us for 
simplicity consider a finite box such that the momenta are 
discrete. (This discreteness in the number of integration 
variables is not related to quantum mechanics - it is the 
same for classical physics.) The integration measure in the 
functional integral ([T]) depends continuously on a. One 
may therefore expect that also the expectation values like 
the average value np(q) depend continuously on a. We 
have just demonstrated that this is not the case for T = 0. 
If for a free theory kp crosses the value q 2 for one of the 
discrete momentum levels the pole in the propagator moves 
from the upper half plane to the lower half plane and np(q) 
jumps by one unit. For stable particles this generalizes to 
the interacting theory. The number of particles with a 
given momentum q is discrete! 

For bosons the situation at T = 0, n ^ is more subtle 
due to the presence of a Bose- Einstein condensate, <f) a ^ 0. 
We now employ the "real basis" of fields \i an d Xi- In 
principle, one has to consider separately the poles of the 
propagator G\\ for the radial mode and of G22 for the 
Goldstone mode. Let us denote by G the 2x2 matrix with 
elements G a b and investigate the inverse matrix P = G _1 . 
By Euclidean time reversal symmetry it has the general 
form 



P 



a, -qoZ^ \ 



(28) 



with a,b,Z^ depending only on q 2 and q 2 ,. This implies 
that g is an even function of qo 

irji\-\ ( b, qoZ^ \ _ 1 a + b 



G={ab + qiZi) 



—qoZ<p, a 
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2 ab + q 2 Z 2 ' 
(29) 

Therefore the location of the poles for Gu and G22 is the 
same, namely for 



(30) 



Here we use that a, 6, Z 2 ^ are real and positive and we as- 
sume Z^ > 0. We recall that a, b, Z^ may depend on q$ 
such eq. (|3H)l is an implicit equation for qo{q)- 

We notice that in the real basis g has poles both in the 
upper and lower half plane. There residua have an opposite 
sign such that both contribute equally to J qo g(q)- If the 
solution of eq. (|30[) is unique we define 

Z (q) = 1 dfl'H^g') = 1 dg^juj,^, 

2 8lO \u (q)>0 2 did \ui (q)«3 

(31) 

This normalization ensures consistency with the complex 
basis in the absence of spontaneous symmetry breaking. 
Relations of the type (fH)|) and (|21[) continue to hold for 
cjq > 0. The residuum of the pole at ujq > is given 



by Z^ 1 , whereas the one for the pole at ojq < has the 
opposite sign. In consequence, both poles contribute an 
equal amount to the qo-integral. This yields for bosons at 
T = and q 2 > 



ZM 



l + 2n{q)' 



(32) 



Thus for single poles we recover eq. (|21l) , while generaliza- 
tions occur if eq. (f3"0")) has more than one solution. 



IV. BOSON OCCUPATION NUMBERS FOR 
LINEAR DYNAMIC TERM 

In this section, we compute the occupation numbers 
for bosons at T = in the approximation that the go- 
dependence of the inverse propagator is linear, as for the 
classical action. This amounts to an ansatz where the co- 
efficients a, b^Zfj, in eq. (|28p depend on q, but not on go- 
For Zj, = 1, a = q 2 /(2M) + 2X^1, b = q 2 /(2M) this 
amounts to the Bogoliubov approximation. We will en- 
counter infrared problems in low dimensions (d — 1). We 
show in the next section how those are cured by a more 
general form of the propagator. In fact, for d = 1,2 the in- 
dependence of the inverse propagator becomes quadratic 
in the regime with nonzero density, where the fluctuations 
of the Goldstone bosons dominate Q. We also encounter 
ultraviolet problems if Z^ differs from one. They are cured 
by a renormalization group improved treatment that we 
discuss in sects. I VII to IXl These issues show that a func- 
tional integral computation of occupation numbers is not 
trivial and needs a thorough understanding of the many 
body dynamics. 

In the limit of go-independent a, b and Z^ one finds 



Z„ = 



2Z 2 \oj \ 



2Zj,V ab 
a + b 



(33) 



In the absence of spontaneous symmetry breaking one has 
a = b and therefore Z u = Z^ (we assume Z^ > 0). In this 
case possible deviations of Z u from one can only arise from 
Z$ ^ 1. In the presence of spontaneous symmetry break- 
ing, however, Z^ deviates from one even for g -independcnt 
a, b and for Z^ — \. We may parameterize 



b = A-, 



2M 



(34) 



where A and may depend on q 2 . This reflects the 
vanishing of a — b for 0o - * and the Goldstone-character 
of Y2, i.e. lim 6 = 0. We will employ later that P is 

related to the second functional derivative of the effective 



action T. The coefficients (|34|) follow from a simple form 
of T, i.e. 



A . 



Z (f> (b*d T cj ) -A<j ) *—cj ) +-X 



They can also arise from a more general form of V. 



(35) 
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With the ansatz J28J fl3j) we find for q 2 > 



I . 1 





(36) 



For = 1 and constant A^ this is the occupation num- 
ber for the Bogoliubov theory. Writing the density n as a 
condensate density </>q and a "particle density" n p 



n p (q) 



(37) 



the Bogoliubov theory yields a finite result for n p only for 
d > 1. Indeed, for small g the occupation number diverges 



q in the Bogoliubov theory 



n p (q) = 



2q ■ 



(38) 



For d = 2, 3 the particle density (|37|) is dominated by 
modes with momenta q 2 ~ AMX^q. 

We will see that quantum fluctuations are responsible 
for important modifications of n p (q) for small q = \q\. The 
dominant effect concerns the go dependence of the propa- 
gator discussed in the next section. In this section we first 
retain the linear dependence of G^ 1 on qo and show that 
the dependence of A^, and Z^ on q does not resolve the 
infrared problems. Evaluating eq. (j3"B"|) for q 2 — ► the 
occupation number diverges according to 



n p (q) = 



ZZfiq) \q]- 



(39) 



The detailed behavior depends on the (f-dependence of A^ 
and Z(j). This can be extracted from a recent functional 
renormalization group study of quantum phase transitions 
For the far infrared behavior q 2 — > the Goldstone 
regime applies. The relevant physics strongly depends on 
the space dimension d. For d — 3 one finds that Z$ vanishes 
logarithmically with q = \q\ 



(40) 



(In the notation of [4| one has Z<f, = S = AS if appropri- 
ate units for length and time are used.) Also vanishes 
logarithmically, A^ ~ Z$. The overall behavior for d = 3 



n(q) 



lnHk 2 /q 2 



(41) 



differs from the free Bose-Einstein distribution which van- 
ishes for q 7^ and T = 0. The infrared contribution to 
the total particle number ~ J dqq 2 n(q) remains finite. 



For d = 1 and d — 2 the issue is more involved, 
convenient to introduce renormalization fields 



0= VA(j) 



It is 



(42) 



The renormalized order parameter fa is related to the su- 
pcrfluid density n s by n s = po = For T = the Galilei 
invariance for non-relativistic bosons yields a Ward iden- 
tity which implies that the total particle density n equals 
Pq. Since A depends on momentum the unrenormalized 
order parameter fa becomes effectively momentum depen- 
dent. In particular, fa vanishes for d — 1 and q — *■ 0, while 
fa remains different from zero. Writing 



AS , S~ q~ vs , A 



(43) 



one finds r/s < — 1 for d = 1,2. For d = 2 the anomalous 
dimension r\ vanishes for q — > 0, while for d = 1 one has 
a range < r\ < 1/2. We also introduce the renormalized 
coupling A = \rj,/A and note that a s = S/X goes to a 
constant for q — > 0. 

The resulting expression for the occupation numbers 
reads for q — > 



n(£) 



1/2 



1 



-1+7)+- 



2^SW)A{q)q 



(44) 



For d = 2 the integral J dqq d ^ 1 n{q) remains infrared finite 
only for r/s > —2, while for d = 1 finiteness requires 77s > 
— 27/. For d = 1 this condition is not obeyed - in this 
case already the "classical approximation" with vanishing 
r\ and rjs leads to a logarithmic divergence. We will see 
in the next section how this too strong infrared increase of 
n(q) is cured by a modification of the go-dependence of the 
propagator. 

We next turn to the behavior for large q as given (for all 
dimensions) by eq. ([56]) . We may use Z<f, = 1 + AZ<j> where 
AZ4, vanishes for q — ► 00. This results in 



n(q) 



M 2 X 2 ^ 



(45) 



q 2 (q 2 + 4MX^ 2 ) 2 

For a first evaluation of AZ^ we may use the flow equation 
of [J], replacing k by jq. For the regime q 2 S> AMX^cj) 2 one 
finds 



d_ 

' dq' 



q—. z 4> 



-(v + vs)z<i> 



32v d 



M 2 X 2 ^{ iq y-* (46) 



with v^ 1 = 2 d+1 K d / 2 T{d/2). In consequence, the leading 
order (with constant X^.fa) 



(47) 



is positive and decreases ~ q . A similar result follows 
from a direct evaluation of the one loop contribution to Z^, 
if Z^ is defined as the coefficient of the term linear in q of 
an expansion of the inverse propagator in powers of qo - cf. 
appendix A. 
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The simple estimate (|45|) for AZ^ cannot be true. The 
decay for large q is too slow. The term ~ AZ^ would 
dominate n(q) for large q and lead to negative occupation 
numbers. Also the g-integral for the total particle density 
would not converge for d > 2. The origin of this problem is 
again the ncglcction of the non-linear go dependence in the 
inverse propagator. It is not justified to evaluate the go- 
dependence by the go-derivative of the inverse propagator 
at qo — 0. If one wants to use an approximation linear 
in g f° r the inverse propagator, should be evaluated 
at the location of the pole rather than near g = 0. We 
will develop in sects. IVII to [X] a renormalization group 
treatment for the computation of n(q). As a result, one 
finds that there is actually no contribution from AZ^ in 
eq. (|45j) . Omitting this term the ultraviolet contribution 
to the total particle number becomes UV-finitc 



The result (|53[) for the low momentum behavior of the 
occupation numbers differs qualitatively from the Bogoli- 
ubov theory ([38]) . It involves the new coupling V, such that 
the coefficient in front of g _1 docs not vanish anymore for 
4>o — ► 0. For g — * one finds [|[ that V — V/A goes to a 



constant such that (cf. eq. 

n(q) 



for d= 1.2 



,-(1-17) 



(54) 



In consequence, the integral J dqq d ~ 1 n(q) remains infrared 
convergent both for d = 2, where r/ = 0, and for d = 
1, where 77 > 0. We note that the infrared convergence 
for d = 1 arises only as a consequence of a nonvanishing 
anomalous dimension 77. The same effect guarantees that 
the correlation function in position space decays with a 
power law for large separation (d — l,r = it) 



M 2 A|$ 



/ f g 2 (g 2 +4A/A ^)' 
as for the Bogoliubov approximation. 



(48) 



V. BOSON OCCUPATION NUMBERS FOR 
GENERAL DYNAMIC TERM 

Let us next address the infrared problem for d = 1. For 
the Bogoliubov approximation the occupation numbers di- 
verge ~ 1/g and the particle density n does not remain 
finite. This problem is exacerbated if the g-dependence of 
\^4>q/Z^ is taken into account. The remedy of this prob- 
lem is related to the go-dependence of a and b in eq. 
that was neglected so far. Indeed, the investigation of 
reveals a crossover to a "relativistic kinetic term" for small 
momenta. For q -C q r one may neglect the off-diagonal 
term ~ qoZ^ in eq. (f2"5)l . Due to the strong decrease of 
for small q (|43|) the go-dependence is now dominated by a 
contribution quadratic in gg in a and b 



V 2M 



2A^0o 



Vq 2 



b = 



2M 



In this region the Goldstone mode dominates 

5 ~ 2 G22 ~ 26 
and the location of the pole is now at 



Aq 2 

2MV 



(49) 



(50) 



(51) 



G ~ (r 2 + 



\ 2 )-. 



(55) 



Away from the infrared limit g — ► one has to follow 
the location of the poles of the general propagator given 
b_y_eqs. (gSJ), They obey, with A = \,p/A , p = 

A(f>l , X Po = X^l, S = Z+/A, 



1 

'W 1 



S z + 2V\ Po + 



±1 {S 2 + 2V\ PQ ) 2 + 2S 2 V q 



M 



(56) 



We discuss here the case where S, V and A are independent 
of go- (A generalization for mildly go-dependent couplings 
will result in implicit equations where these quantities have 
to be evaluated on the location of the poles.) On the other 
hand, the couplings depend on q. For large q the fluctua- 
tion effects are small, such that V — * 0, S — > 1, while for 
small q one has 5* — > 0, V 7^ 0. All solutions of eq. ([5^)1 
correspond to negative q$ and therefore to pairs of purely 
imaginary go. 

Let us start with V = where the poles at 



g = ±iS~ 




(57) 



With 



can be associated with the Goldstone mode. We observe 
positive and negative frequencies for the excitations, cor- 
responding to particles and antiparticles, the latter being 
associated to "holes" in the condensate. For small nonzero 
V, the poles (|57|) get slightly shifted. Two new poles ap- 
pear at large |go|, 



[2AV\ 1/2 
Z^ = 2^M - y — ) q 

the particle density for g — > (g 7^ 0) obeys 



(52) 



g = ±i 



V 



(58) 



As S/V decreases towards zero, these modes correspond to 
the radial modes of the relativistic model, i.e. for S — * 
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The Goldstone poles take for S — > the standard relativis- 
tic values 



9o,G 



q 2 



2MV 



1 



S 2 



(60) 



The radial mode corresponds to collective fluctuations re- 
lated to the size of the condensate po, whereas the Gold- 
stone mode is associated to phase changes of the conden- 
sate. More generally, we refer to the negative (positive) 
sign of the square root in eq. ([56]) as the Goldstone (ra- 
dial) poles. 

For a computation of the occupation number with the 
propagator P5|) we evaluate for q ^ 



h(q) 



qo J qo 



1 



2j qo ab + Z 2 q 2 

i f ( — 



\2M 



+ Vq% + \p )det- 1 {q) (61) 



with 



det(g) = 



\2M 



Vq 2 + 2Ap 



V 2M 



V 2 U + JL( a + y/0) 



(62) 



and 



a = S 

(3 



2VXp + 



Vq 2 
M 



(S 2 + 2VX Po ) 2 + 2S 2 Vjj. 



Writing 

det _1 (g) 



V M \ % + IV' 1 % + 2V i ■ 



(63) 



(64) 



we recognize the contribution from the radial and Gold- 
stone poles 



h(q) 



1 



J 9o \9 + 2V 2 



<?0 2 



2V a 



(65) 



The ^-integration is performed easily 



h(q) 



1 



2^^ 
M 
2V2AV 



i-ff 2 /V7? | i + ff 2 /v/7 



{g 2 (g 2 + 4MAp )} 



-1/2 



(66) 



with 



T 



a 



(67) 



For the limiting case V — > one has a = y 7 /? — > >S 2 + 
V% 2 /A/ + 2Ap ),^ -» V2(V/S)(q 2 /M + 2X Po ) and we 
recover eq. (|36[) . In the limit S* — > we can use a = 
(q 2 + 2M\p Q )V/M, y/]3 = 2V\p and T = ^/V/M(^f + 

q 2 + AMXpo) such that we recover ([53")) . 

For an estimate of the leading behavior for q — > 00 we 
need J-/V and therefore the g-dependence of V. A renor- 
malization group estimate Q yields 



lim V(q) 

q — >oo 



d(6 - d) 



AW 



d-6 



(68) 



In consequence, the terms ~ Vq 2 in eq. ([6"3"| become sub- 
leading and the ultraviolet behavior is well described by 
the linear dynamic term in sect. HVl 

In the spirit of a renormalization group improvement, 
which will be motivated in the next sections, we may derive 
a differential equation for the g-dependence of n p 



dn p (q) _ dh(q) 



dhiq 



dh\q 
M 



2V2AV 



7(9), 



{q 2 (q 2 + 4M\p )} 



-1/2 



q 2 + AMXpo 



T 



dT 
d In q 



(69) 



Here dJ-/dlnq should be evaluated for fixed couplings. 
The advantage of such an equation is that the ultraviolet 
problem associated with the additive counterterm is now 
absent. For large q the function 7(g) decays ~ q~ 4 . We can 



therefore integrate eq. ([69]) . starting with an initial value 
for large q = A for which the perturbative result eq. ([3"6"[) 
with Zff, = 1) can be trusted. Furthermore, the "renormal- 
ization group improvement" of using g-dependent couplings 
A, S, A^, (/>oj V will make the solution of eq. ([6"9"[) rather ac- 
curate. It can be used for all dimensions and is free of 
infrared or ultraviolet problems. 

Eq. ([nU) can directly be used to compute the particle 
density or condensate depletion n p = J^n p (q) by solving 

a second order differential equation. With x = ln(q/fc) 
and k a suitable momentum unit, one may compute n p = 
fe n (a; — > —00) from the differential equation 



d 2 " ,9n dx , . 



(70) 



using the initial condition n(oo) = (or a similar condition 
for x = ln(A/fc)). 



VI. RENORMALIZATION GROUP 
IMPROVEMENT 

In view of the rather complicated momentum depen- 
dence of the full propagator it is useful to consider the 
contributions from fluctuations with different momenta p 
separately. This will also allow us to use a renormaliza- 
tion group improvement by employing effective couplings 
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appropriate for a given momentum range. In this section 
we will follow a heuristic approach. The second part of 
this paper will present a more formal development, based 
on exact functional renormalization group equations. 

We start by adding to the inverse propagator an infrared 
cutoff piece such that 



G k (q ,q) = (P(qo,q) + Rk(q)y 



(71) 



The infrared cutoff is a 2 x 2 diagonal matrix and suppresses 
all modes with q 2 < fc 2 . We choose 

Rk(q) = ^(k 2 -q 2 )9(k 2 -q 2 ) (72) 

such that the modes with q 2 > fc 2 are not affected. On the 
other hand, the increase of the classical free propagator for 
small momenta, G ~ q~ 2 , is stopped at q~ 2 = k~ 2 . We 
define a fc-dependent occupation number (for q ^ 0) 



n P*@ = \J TrG k {q)- 1 - 



and consider the derivative 



dkn Pik (q) 



Tr{d k R k (q)G 2 (q ,q)}. (74) 



We note that d k Rk vanishes for q 2 > fc 2 , such that the 
running stops once k < q. Furthermore, for the classical 
propagator the r.h.s. of eq. (|74|) becomes independent 
of q for q 2 < k 2 . The go-integration on the r.h.s. of 
eq. (|T4[) will be dominated by fc-dependent poles. The 
physical occupation number obtains by solving eq. (174|) 
for k — > 0. Indeed, the absence of the cutoff for k = 
implies n p (q) = n pfl (q). 

More explicitely, the flow equation reads 



dkn Ptk (q) 



kA 
2M 



6{k 2 -q 2 



»(-!) 



{G\i + G\ 2 + 2G12G21) 



We will neglect the anomalous dimension 



V 



~ k Wk lnA 



(75) 



(76) 



in this section. Inserting the ansatz (|28|) and approximat- 
ing a, b and Z^ to be go-independent one finds 



scale fc, such that they approximate best the relevant cou- 
plings in the momentum range p 2 rs fc 2 . We will moti- 
vate this procedure more explicitely in the second part of 
the paper. A crucial property for the ultraviolet issue is 
Z<f,(k — > 00) = 1. This simply follows from the absence 
of fluctuation effects for an infinitely large cutoff. This 
property guarantees the correct initial value for the flow 
np.k^oo(q) — for all q. Indeed, the classical contribution 
to n v _k for very large fc yields eq. (|36|) with q 2 replaced by 
fc 2 and all parameters evaluated at the scale fc. It is obvious 
that Z${k — > 00) = 1 is needed for a vanishing occupation 
n Pi fc in this limit. The "initial value" Z$(k — > 00) = 1 is 
well compatible with a nonzero value Z^,(k — > 0). While 
only the latter is used in eq. ([56"]) . the flow equation ([71)1 
involves the value of Z^(k) for all fc 2 > q 2 . We will see 
that the approach with a varying Z^(k) avoids in a simple 
way the difficult ultraviolet issues encountered in the direct 
evaluation of eq. ([56]) . 

Let us consider the range of high momenta q 2 3> 
AMX^q. In the relevant fc-range fc 2 > q 2 we can ap- 



(73) proximate 



d k n p (q) = -4M 2 A 2 (fc)^(fc)Z- 1 (fc)fc~ 5 0(fc 2 - q 2 ). (78) 
This yields 

n p (q) = Af 2 (A 2 ^z- 1 )g- 4 (79) 
where we have defined the weighted average 



(A 2 ^ 4 ^ 1 ) = j dk\ 2 (k)4>Uk)z^(k)k- 5 1 J dkk~ 5 . 

g 1 

(80) 

The factor fc 5 strongly enhances the weight of the range 
fc q and we may approximate 



(81) 



Comparison with eq. (f3"6"|) . ([43]) yields AZ^ = and there- 
fore effectively Z$ = 1 for the range q 2 3> 2MA 2 ,</> 2 ,. No 
difficult ultraviolet issue needs to be addressed anymore. In 
particular, the renormalization group improvement avoids 
the need to evaluate the propagator at nonzero qo, corre- 
sponding to the location of the pole. We have also found 
that the parameters \$, Z$ appearing in eq. (j3"B"|) should 
be replaced by g-dependent effective couplings. We empha- 
size that our argument in favor of AZ<p(q) = is valid for 
q 2 3> 2MAj,</> 2 and does not necessarily apply to small q 2 . 



kd k n Pi k(q) = - 



4M 2 A?^ 

-6{k 2 -q 2 ). (77) 



Z^kik 2 + 4MX^ 2 )V 2 



VII. OCCUPATION NUMBERS FROM SOURCE 
TERMS 



Using for fc — > 00 the initial condition n p .fc^oo = an 
integration for fc- independent X^,,<po,Z^ yields for fc = 
eq. ©■ 

The renormalization group improvement uses fc- 
dependent parameters A^, 0o and Z^. These fc-dependent 
couplings are in turn computed in presence of the cutoff 



In sects. IVIIUlXl we address the formulation of the flow 
equations for fc-dependent occupation numbers more for- 
mally. We will derive an exact functional differential equa- 
tion and propose non-perturbative approximations. The 
solution of the flow equations offers an alternative to the 
direct computation of n(q) from the propagator (|T2")) . The 
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go-integration in eq. (|12p depends in a critical way on the 
precise knowledge of the go-dependence of the propagator 
which is quite involved. In contrast, the flow equations are 
less sensitive to this dependence. Only a small go range 
influences the running at a given renormalization scale and 
rather crude approximations are often sufficient. The prize 
to pay is the numerical integration of the flow equations. In 
this paper we will not perform detailed quantitative com- 
putations. We rather use the flow equations for a justifica- 
tion of the renormalization group improvement discussed 
in sect. [VTl 

In the present short section we set up the formal tools 
for the computation of the momentum distribution n(q). 
For this purpose we add appropriate source terms to the 
action |T]) 



s ->• s - h(q)x*(q)x(q)- 

Jq 



(82) 



The partition function Z depends then on h{q) and we can 
formally extract n(q) as 



(83) 



Here we set the source term to zero after the differentiation 
and the minus sign applies for fermions. We note that the 
source multiplies a composite operator and is nonlocal in 
the fields x( x ) 



h{q)x* '(q)x(q) 



h(x- y)x*(T,x)x{r,y), 



(84) 



t J x J y 



with h(z) the Fourier transform of h(q). Eq. I|84p shows 
that the source term remains translation invariant. 

We will conveniently work with the effective action T[(f>] 
which generates the IP I correlation functions. It ob- 
tains by introducing local linear sources j(x) for x( x ) an d 
performing a Legendre transform of lnZ[j], with cf>(x) = 
(x( x ))\j m the presence of sources. The nonlocal sources 
h(q) are not affected by this transformation and T depends 
now on h(q) and 4>(x). We also add a field-independent 
piece linear in h(q) which takes care of the "additive renor- 
malization" ±1/2 on the l.h.s.of eq. 



T[$\ = -\nZ[j] 



(85) 



h(q). 



lq J q 

It is easy to verify that the occupation numbers obey 

1 ST 



n(q) = - 



Q d+1 Sh{q)^ o{x)h=0 



(86) 



Here the partial derivative S/5h(q) is taken at fixed <j)(x) 
and evaluated for the solution of the field equation <pg (x) 
which obeys 



sr 



o. 



(87) 



Our goal will be the computation of T[4>; h] in linear order 
in h. Expanding 



r[0; h] = T[4>] - / Affa q\h(q) + 0(h 2 ) 

Jq 

the occupation numbers follow directly from J\f 



n(q) 



1 



d+l 



(89) 



The field equation for </>q only involves r [<f>]. 



VIII. EXACT FLOW EQUATION FOR 
OCCUPATION NUMBERS 



The difference of Af[4>; q\ from its classical value results 
from quantum and thermal fluctuations. We include these 
fluctuation effects stepwise by introducing first an infrared 
cutoff R k which suppresses the fluctuations with momenta 
q 2 < k 2 . (For fermions we may choose to suppress the 
fluctuations with \q 2 — kp\ < k 2 .) At the end we lower k 
to zero such that all fluctuations are included. We achieve 
this strategy by adding to the action (fT]) an infrared cutoff 
term Q. 



Rk(q)x*(q)x(q)- 



(90) 



In turn, the effective action is now replaced by the average 
action Tfc which depends on k [||, Q- With Rk(q) diverging 
for k — > oo all fluctuations are suppressed in this limit and 
one has Tk^oo = S. On the other hand Rk{q) = for 
k —* implies T^o = T. The average action therefore 
interpolates smoothly between the classical action for k — > 
oo (or k = A) and the effective action for k — ► 0. 

Its dependence on k obeys an exact flow equation @ 



1 



Tr{d k K(T {2) [0; h]+TZ)- 1 }, (91) 



with H(q,q') = Rk(q)S(q — q'). The second functional 

(2) 

derivative T k is given by the full inverse propagator in 
the presence of "background fields" </>. For a homogeneous 
background field one has 1^ (q,q r ) — P(q)5(q ~ q), with 
P a matrix in the space of complex fields (0, </>*) or real 
fields (0i, 02), with, <f> = -^(4>i + i<t>2)- The trace involves 
a momentum integration and trace over internal indices. 
We will choose here the basis (<f>i,(f>2)- For homogeneous 
background fields (f> eq. ([91) takes the explicit form (with 
tr the internal trace) 



d k T k 



d+l 



tr / d k R k {q){ l2 {q) + Rkiq)' 1 ■ (92) 
Jq 



The precise shape of the cutoff function R k is, in principle, 
arbitrary. As an example, one may consider the "optimized 
cutoff" [13] 



' 2 )&(k 2 



q 2 )- 



(93) 
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We observe that for Rk depending only on q the nonlocal 
source term ([H2"]l may be viewed as a shift R k (q) — ► Rk(q) — 
h(q) such that a variation with respect to h can also be 
interpreted as a variation in the space of cutoff functions. 

The presence of the nonlocal source h(q) does not change 
the general structure of the exact flow equation (f9Tj) . Em- 
ploying again the expansion (|88|) (now with all quantities 
depending on fc) the flow for T k [4>] obtains from eq. (|9"Tj) by 
simply omitting the /i-dependence, corresponding an eval- 
uation for h — 0. The solution of the field equation ([57]) 
4>o(k) will now depend on k. The exact flow equation for 
Af k [4>] q\ obtains from eq. (|9ip by taking a derivative with 
respect to h(q) 

dkJVk = ~Tr{d k H(T^ +K)- 1 Mj i 2 \T£ ) +K)- 1 }. (94) 
Here Afjf' (q' , q") is the second functional derivative 

<^ (95) 

Eq. (|94p will be our basic starting point for a computation 
of n(q) through a flow equation. We will consider a fc- 
dependent n k {q) as defined by eq. (|89[) and extract the 
occupation numbers for k — > 0. 

For a homogenous (po(k) one has A/J (q',q") = 
V2(q')&(q' — q") and the flow equation for Afk simplifies 

dk^So] = (96) 
-^T tr j ql d ^k (q >2 (q')lP(q') + Rk (q')} ~ 2 

where v 2 is evaluated for = <j) . This results in the final 
flow equation for the occupation numbers 

dknk(q) = - \ tr J dkRk(q")v2(q']q) 

[P(q>) + R k{ q>)]-* + ?10-d k %. (97) 

Here the second term accounts for the fc-dependence of the 
solution <po(k). 

For k — > oo the fluctuation effects can be computed ex- 
actly. For bosons and a < we start for k = A with the 
initial value 

and 

A/l 2) [Ml (</,</') = ?'W-<f')Sab. (99) 

More precisely, the initial value for k — * oo has a classical 
term and a one loop contribution. For bosons and er < 
the latter precisely cancels the classical contribution from 
the last term in eq. (|55|) . Since there is essentially no 
running between fc->oo and k = A we may begin the flow 
at k = A. (For fermions and a > the one loop effect is 



less simple, see eq. (|140l) below.) According to eq. (|98|) 
one has for a < the initial value 

n A (§) = 4> 2 (A)5(q). (100) 

The only possible contribution arises from a condensate 
4>q (A) ^ 0. For g^Oa nonzero occupation number has to 
be built up as a result of the flow. 

We want to extract n(q) from a solution of the flow equa- 
tion (|94|) for k — > 0. This will only be possible in suitable 
approximations. One may imagine to solve first the flow 
equation for T k and then use the result in eq. (|95[) . For 
this purpose we will need ^2(9'; <tj- 

The quantity v obeys a type of Ward identity. We have 
defined n p (q) in terms of the fc-dependent propagator. Us- 
ing the exact flow equation for the exact propagator Q 
yields a second exact equation for the flow of n p . The com- 
parison of the two exact flow equations yields the 'Ward 
identity" . 

Let is demonstrate this for a homogeneous setting where 

If- 1 

n P ,k(q) = r/ trG(q ,q)--, (101) 
1 Jqo 1 

with G the fc-dependent propagator matrix. For q ^ this 
yields the flow equation 

d k n k (q) = Uq) + ^#^4 (102) 

with 

Uq) = \l trd k G(q ,q) (103) 
= ~\ I tr[(P + R k )- 2 dk(P + Rk)](qo,q)- 

1 Jqo 

If we decompose 

u(q'; q) = 5{q' - q) + u v (q' , q) (104) 
we obtain by comparing eqs. (|9"7| and (|102[) 

tr j d k R k {q')v p {q'-q)G 2 ^)^tr f d k P(q , q)G 2 (q , q). 
Jqi Jq 

(105) 

With F^ 2 '(q, q') — P(q)6(q — q') the exact flow equation for 
P follows from the second functional derivative of eq. (f9"Tj) 
6]. It can be written as a one loop expression 

d k P(q)=tr[ d k Rk(q , )G 2 (q , )H(q , ;q), (106) 
Jq' 

with Ti a matrix both in the space where G acts (and over 
which the trace in eq. (|106|) is taken) as well as in the space 
where P acts. We identify 

v v {fl\ q) = tr f H(q'; q Q , q)G 2 (q , q), (107) 
Jqo 

with tr denoting the trace in the space where G acts. 
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IX. QUADRATIC TRUNCATION 

The flow equation ([97]) is exact. However, even if P(q') 
would be known from a solution to the flow equation for 
r^.[0], the r.h.s. still involves the unknown second func- 
tional derivative One may derive an exact flow equa- 
tion for V2, but in turn it will involve even higher functional 
derivatives of AT. We need to approximate the solution by 
a suitable ansatz (truncation) for the general form of N ' ■ 

As a first approach we consider bosons and truncate AT 
in the number of fields. We keep only terms with up to two 
powers of <f> 

1 a J l'o 
Jq' 

-n d+1 Z Ptk (0,q)4>l (108) 

The initial values are 

n P ,K{q) = , Z CtA (q' ,q) = 1 , Z PtA (q',q) = 0. (109) 

For 4>o(q) = (f>od(q') the momentum distribution reads 

n(q)=n p ,o(q)+Z c , (0,0)4>l6(q). (110) 

The condensate density 

n c = Z cfi (0,0)4> 2 (111) 

differs from the bare order parameter 0q only if Z c (0, 0) ^ 
1. The total density n and condensate fraction f2 c read 



n p (q) , VL C = n c /n. 



(112) 



We need to understand the meaning of the functions Z ck 
and Z Pt k and start by defining them properly. Let us first 
use a combined "density transfer function" Z^ 

Z h (q', q) = Z p (q', q) + Z c (q' , q)8{q " - q) (113) 

and define Z p and Z c by their behavior for q ' — > q. We 
want to choose Z p (q' , q) such that it remains finite for q ' = 
q. We therefor define 



Z c(%,q) = ^ d 1 J Z h (q ,q)8{q' - q), ( 



114) 



which yields Z p implicitely by subtraction in eq. (|113[) . 
(One may replace S(q ' — q) by a smoothened ^-function 
and change Z p {q',q) correspondingly. This concerns the 
practical definition which particles are counted as having 
zero momentum.) 

In terms of Zh the second functional derivative reads 

A^ 2 > (<?',<?") = Z h {q\q)5{q' - q") , v 2 {q> \q) = Z h (q' ,q). 

(115) 



All quantities need to be evaluated for the solution of the 
field equation <po(k). The substraction of a constant term 
~ Z Pi k4>o m the truncation (|108p guarantees that the ex- 
pansion in powers of fields is correctly performed around 
4>o(k) and that Z p does not contribute to the particle den- 
sity. The flow equation for n{q) ([9"T]) has both a contribu- 
tion from Z p and Z c 

d k n{q) = --tr [ d k R k (q)Z c (q' ,q) 

[P(q^q) + R k {q)Y 2 
-\tr £ d k R k (q')Z p (q>,<f)[P(q>) + R k (q')]- 2 



(116) 



In the remainder of this section we argue that the con- 
densate density is precisely given by (f> 2 . This means that 
Z c (0, 0) = 1 is not renormalized, as will be shown more 
precisely in the next section. We observe that both con- 
tributions to eq. (|116p from Z p and Z c are regular for 
q 2 — > and therefore do not contribute to the running 
of the condensate density n c . The latter simply obeys for 
0o ^ 



(117) 



and d k n c = for 0o = 0. 

We will still need the (^-dependence of n(q). Again, we 
derive a flow equation for 



dn k {q) 



It reads = +K) 

d k a(q) = -^—Tr {d k K^{QN^Q) 



1 

-tr 
2 



d+l 

Gd k R k G 



(118) 



(119) 



- GZ h +Z h G - 



where we have omitted a term ~ dZ^/d^ in the last equa- 
tion. The only dependence on q arises from Zh and we 
observe that the r.h.s. of eq. (] 1 1 9[) remains regular for 
q 2 — > 0. In consequence, we can use for all k a non-singular 
&k(q) 

a(q) = 5{q) + a{q) (120) 
and the condensate density obeys 

= 1 , d k n c = d k 4>l , n Cik = <f>l(k). (121) 

O(f> 

This closes the argument that the condensate density is 
not renormalized and shows that n p (q) has no contribution 
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~ S(q). For the SSB regime with <fio(k) 2 > we may define 
Z c (0, 0) = n c /4>l such that 



a fe z c (o,o) 



dZ c (0,0) , 2 



(122) 



A simple solution with the initial condition of a <fr indepen- 
dent Z C)A (0,0) = 1 is Z c (0,0) = 1 for all k. 



X. DENSITY TRANSFER FUNCTION 



The right hand side of the density flow equation ()116[) 
involves the "density transfer function" Zh(q',q) (|113|) . In 
order to understand its role we neglect for a moment the 
dependence of Zh on q' . The first two terms in eq. (|116p 
can then be written as 



8 k n{q) = Z h {q',q)d k n{q ') + ... (123) 

Jq ' 

where dtfi(q) accounts for the contribution of a loop with 
momentum q' 



d k h(q') = -\tr f d k R k (q')(P(q ,q') + R k (q')) \ 



(124) 



For a go-dependent Z% the folding (|123[) is in a certain sense 
averaged over q' . The density transfer function describes 
how a loop with momentum q ' influences the occupation 
number n{q). 

The flow equation for Zh{q' ,q) obtains by taking the 
second functional derivative of eq. (|M|) with respect to (j). 

In the truncation (I108P Af k does not depend on <j) and one 

finds with g = (r£° +TZ)- 1 

d k fV^(q',q") = l -Tr{[gd k KQN^Q 
^ r (2) 



+gu {2) gd k ng] 



5<l>l(qi)84> b (q») 
6TW 



-\TAgd k ng^—gN^g^ 

2 T Wb(q") Scj)* a (q>) 

-i ^g(N {2) gd k n + d k ngN {2) )g 

dT (2) dT (2) dT (2) Q T (2) 

-g — , „ N + „ T ; ,,\ g r\~ .,. r » ) \ • (125) 



KW) d4> b { q ") dUi") dKW) 



The first term involves a quartic vertex and the two last 
terms arise from cubic vertices present for ifio ^ 0. We will 



approximate the vertices here by pointlike vertices 



(2) 
cd 



*r£V,p") 



Hb{q") 

For a quartic coupling one has 



= KbcdS{p' - p" + q' - q" 
= lacd<j>o5{p' ~P" + ?')) 

= TbcdfaStf -p" -q"). 



(126) 



lacd = lacd = V2X{S a l5cd + S c lS a d + SdlS a c), 
Xabcd = ^(SabScd + 8 ac Sbd + Sad^bc)- (127) 

For a homogeneous setting g is diagonal in momentum 
space. One obtains the flow equations 

2 r 

d k N£ (<?', q"; q) = 8{q' - q") £ / d k R k (jf) 

a.b=li, 



{\ccab(G 3 )ba(p')Z h (p , q) 

-f E {(G 2 u( P ') 



d,e=l 



(lcbdJcea(G 2 )de(p' + q')Z h {p' + q' , q) 

+lcbdlcea{G 2 )de{p' - q')Z h {p' - </, q)) 

+ 2(G 3 ) ab {p')Z h (p', q){ lcbd lceaG de {p' + q') 

+lcbdlceaG de {p' ~ g'))]} (128) 

Since d k M^ gets, in principle, also contributions from 
higher order terms neglected in our truncation we identify 

d k Z h (q', q)5(q' - q") = d k Afg> (q' , q"; q). (129) 



The r.h.s. of eq. (|128[) remains finite for q' — > q. The 
contribution from Z c changes the g'-integration to an inte- 
gration over q' with values of p fixed by the 5- function in 
Zh, i.e. pi = q or $ = <f±g '. The go-integration is well be- 
haved due to the infrared cutoff R k in G~ 1 . This leads to 
the important conclusion that renormalization effects for 
the condensate contribution are absent 



d k Z c (q / o ,q) = 0. 



(130) 



The condensate density is therefore given by the "bare" 
order parameter 



n c = 4> t 



(131) 



This is consistent with the discussion after eq. (|122[) . 

In contrast, the build up of the occupation number for 
uncondensed particles is influenced by the renormalization 
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explicite cutoff (|93|) such that 



d k n p {q) = 

- \tr [ d k R k (q) [P(q' : q)+R k (q)}- 2 
1 J l'o 

™*r I d k R k (q r )Z p (q',q)[P(q') + R k (q r )}- 



kd k n(q) 



9(k 2 ~q 2 ) 
k 2 -q 2 



-2 

o \ (136) 



The occupation number n(q) flows only as long as A: 2 > q 2 , 
then the flow stops. In the range k 2 > q 2 the flow equation 



(132) simplifies 



Only the first term in eq. (|132p is present in the "classical 
contribution" that neglects the renormalization effects for 
Zh- For the second term the initial value Z p = does not 
remain a solution of the flow equation. A nonzero Z p is 
generated by the contribution <5((f' — q) in Zh(q' , q) in the 
r.h.s of eq. (|128[) . Whereas the first term involves only "on 
shell momenta" where the propagator in the loop is fixed 
to q ' = q , the second term gets also contributions from "off 
shell momenta", i.e. from q' different from q. This reflects 
the broadening of the density transfer function Zh(q',q) 
due to the renormalization flow. 



XI. OCCUPATION NUMBER FLOW AT LOW 
TEMPERATURE 

In order to get some practice with the flow equations we 
first check simple situations. We start with the vacuum. 
The term ~ 0q in eq. (| 1 28(1 is absent. We want to compute 
the contribution from Z c to the flow of Z, D 



d 

dk 2 



1 

2M 



2M 



= 



(137) 



Since for large k = A the initial value n p ^(q) vanishes 
we conclude that only the codensate n c contributes to the 
particle number, n = n c . 

For fermions a positive chemical potential a > is rel- 
evant. Now eq. (|136[) is ill defined when max (q 2 ,k 2 ) = 
2Ma. This reflects that the cutoff (j9"3"|) is not well adapted 
to the problem. Indeed, instead of cutting out the small 
momentum fluctuations an efficient cutoff should rather 
suppress the fluctuations close to the Fermi surface, i.e. 
with \q 2 — 2Ma\ < k 2 . For the fermions we will therefore 
employ the cutoff 



R 



(F) 



= ^(k 2 sgnZ-O0(k 2 -O , 



[a a. 



(138) 



(F) 

For positive a the propagator in presence of R k becomes 
in the vicinity of the Fermi surface 



2 r 

d k Z p (p', q)=Y^ *22aad k R k {q) / (G 3 ) aa (^, q). (133) 

a=l Jlo 

It is convenient to work in the complex basis with diagonal 
G. For large |g | the integrand decays ~ |<Zol~ 3 such that 
we can close the integral either on the upper or lower half 
plane. The integral vanishes, as most easily seen in the 
complex basis where G has only one pole in the upper half 
plane and we can evaluate the closed integral in the lower 
half plane. In consequence, Z p stays zero for all k and the 
density transfer function remains a <5-function 



Z h {p' = 5{jf - q). 



(134) 



Inserting this result in eq. 
tributes 



(|132|) only the first term con- 



d k n p (q) = -\d k R k {q)tr ( G 2 (q' ,q). 
2 J i' 



(135) 



and 



Again, G {q^q) decays sufficiently fast for large |<7q| 
the r.h.s. of eq. (|135p vanishes, as it should be. 

Next we may have a look on the free theory for T = 0. 
In the absence of interactions the density transfer func- 
tion does not flow and eq. (|135[) remains valid, with 
the classical inverse propagator. Let us use for a < the 



G 



k 2 



( F) 

such that the pole in presence of R k remains in the same 



(139) 



half plane in complex go-space as for R 



0. As a result, 



the initial value of the flow is given by the Fermi distribu- 
tion 



n F Aq) = 0(2Ma-q 2 ). 
In turn, the flow equation (|137|) is replaced by 



(140) 



d 

dk 2 



n F {q) 



sgn(q 2 - 2Mo) 
2M 



'/.) 



J 9o 



2M 



sgn(q 2 - 2Ma) 



(141) 



and the Fermi distribution is maintained for k — > 0. For 
T > the r.h.s. does not vanish and the flow accounts for 
the smoothening of the Fermi distribution. 

We next turn on the interactions. The initial values 
n p,h{q) = , nF,\(q) = 0(2Ma — q 2 ) are not changed. 
However, the flow of n(q) can now account for nonvanish- 
ing n p (q) or for a smoothening of the Fermi surface due to 
fluctuations. There are two possible effects. The first con- 
cerns the modifications of the propagator. Furthermore, a 
second effect arises from the build up of a nonvanishing Z p . 
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Let us first consider bosons and employ the cutoff (|93 
The flow of n(q) (cf. eq. Q124p becomes (77 = — kd k \nA) 



d k n(q) 



kA 
2M 



6(k 2 



'>('-?) 



(G11 + ^22 + 2G12G21). 



(142) 



% 



We use the general form 



Gu + G22 + 2G12G21 — 



(a + R k ) 2 + (b + R k ) 2 -2g 2 Z 2 
[(a + R k )(b + R k ) + q 2 Z 2 " 



(143) 

with momentum dependent coefficients a, 6, Z^. In order to 
understand the general structure we concentrate on large 
q 2 where perturbation theory should remain valid. We 
want to see in which limit we recover the Bogoliubov theory 
and what are the extensions of it. For this purpose we may 
and employ the ansatz (f34|) . Neglecting rj and a possible 
go-dependence of A, A^, Z^ this yields [Z^ = AS) 



(k 2 /2M + 2\^ 2 ) 2 + (k 2 /2M) 2 - 2q 2 S 2 



[(fc 2 /2Af)(fc 2 /2Af + 2X^ 2 ) + q 2 S 2 Y 

-8(k 2 -q 2 ). (144) 



4M 2 A 2 ^ 



Z^k 2 (k 2 +AM\^ 2 y/ 2 

For fc-independent X^,4>q and Z$ the integration of eq. 
(fTl4"|) yields 



h(q) 



dk 



4M 2 A 2 ^ 



Z fc 2 (fc 2 +4MA^ 2 ) 3 / 2 



1 

2Z, 



1 + 



4M 2 A 2 «^ 



h A (q) (145) 
1/2 

I -1 



q 2 (q 2 +4MX^ 2 ) / 

-{q^ A} + h A (q) 

For A — > 00, n A (q) — (we only consider here q 7^ 0) the 
sum of the last two terms vanishes, such that 



h{q) = 



2Z d 



4M 2 \ 2 $ 



1/2 



q 2 {q 2 + AM\^ 2 ) 



- 1 



(146) 



If we replace h by n this is almost eq. (f31)f - the differ- 
ence vanishes for Z$ = 1. The Bogoliubov result there- 
fore obtains if we use the "classical transfer function" 
Z c = 1, Z p = and neglect the dependence of a,b, Z<f, 
on qo,q and fc. 

Extensions beyond the Bogoliubov theory take into ac- 
count the running and the momentum dependence of the 
parameters as well as a non-trivial transfer function. In 
a first step we may use in eq. (|146|) effective parameters 
X^,(j)o, Zc/, which depend on q 2 . Based on an approximate 
solution of the flow equation for the occupation numbers 



we argue that one should use the values X^iq, k 2 — q 2 ) and 
4>o(k 2 = q 2 ) instead of the values at k = 0. This is justified 
since the build up of a nonvanishing n(q) is dominated by 
the range k 2 ~ q 2 . 

Indeed, the flow (|144|) stops for k 2 < q 2 , while for 
k 2 3> q 2 it is suppressed by inverse powers of k. Using 
also Z<f, = Z,p(q, k 2 — q 2 ) we note that Z§ multiplies now a 
term that vanishes for q — > 00, in distinction from eq. (|36|) . 
As already noted in sect. IVII this renormalization group 
improvement gets rid of the dangerous term ~ AZ^ in eq. 
(|4"5)) . Indeed, the replacement of Z ( f > (k) by a /c-independent 
constant is only allowed if Z<p multiplies a quantity that 
vanishes for large q, as in eq. (|146|) . This is necessary since 
Z^ multiplies in the flow equation Q144p a term ~ k~ 5 for 
large k. The properly implemented renormalization group 
flow takes care of the ultraviolet problems, even if the reso- 
lution of the go-dependence of the propagator is only poor. 
This has direct consequences for the A-dependence of the 
relation between particle density and chemical potential, 
n(cr). For large A the momentum integration (for q < A) 
of n(q) yields n ~ A d ~ 4 , with a coefficient that depends on 
a. For T = the cutoff dependence of n(a) is suppressed 
by (2Mcr/A 2 )( 4 - d )/ 2 . 

A second improvement beyond Bogoliubov theory be- 
comes necessary if we consider small q. Now the flow of the 
occupation numbers should also contain the term ~ Vq 2 in 
the inverse propagator, as we have discussed in sect. HVl 

Let us finally discuss the flow of the density transfer 
function which generates a nonvanishing Z p according to 
eq. (|128p . The cubic couplings must have an even number 
of indices taking the value two. They are given by 

7212 = 7221 = 7122 = 7212 = 7221 = 7122 = V^A^,, 

7111 = 7111 = 3\/2A . (147) 
and the flow induced by these cubic couplings becomes 

d k z p ^,q ) = -^fj pi e(k 2 -p 2 ) 

{ [{G 2 )n{p'){G 2 ) 22 {p' + q') + (G 2 h 2 (p'){G 2 )n{p' + q') 

+G 2 i 2 {p')(G 2 )i 2 {p' + q') + {G 2 ) 2 i{p')(G 2 ) 2 i{p' + q')] 

Z h (p' + q',q) 
+ [(G 3 )n(p')G 22 (p / + q') + (G 3 ) 22 (p r )Gn(p' + q') 

+ (G 3 )i 2 (p')Gi 2 (p' + q') + (G 3 ) 2 i(p')G 2 i(p' + q')] 

Z h (p',q) + (q' ^ -<?')}• (148) 



We evaluate eq. (|148l) for q' = 0, with the ansatz ((28)) and 
a = a + R k , b = b + R k . 



d k Z Pl3 (0,q) = -4A 2 fc0 2 / 0(k 2 -p 2 )Z h (p',q) 

Jp' 

(H + b) 2 (ab + p> 2 Z 2 )- 3 . (149) 

For the high momentum region we use eq. (|34[) and the 
cutoff R k 172]) such that 



A 

b =2M k2 . a = b + 2A\ 4> 4> 2 . 



(150) 
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The flow 



d k Z p , 3 (0,q) 



12M 2 X 



2\2 ±1 



%<$>t (k 2 + 2MX^ 2 ) 2 



Z^A 2 ^ (k 2 + 4MA^g)5/2 
[o{k 2 -q 2 ) + [8(k 2 ~f)Z p (p 7 q)} (151) 



stops for q 2 <k 2 and is dominated by A: 2 
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As long 



as Z p remains small we can neglect the second term in the 
curled bracket. Omitting also the contribution from the 
four point vertex one finds for large q 



Z P (0,q) = 



3M 2 X 2 J 2 



(152) 



We conclude that Z p appears only in higher order in per- 
turbation theory (~ A^) and is suppressed for large q. For 
a discussion of the ultraviolet behavior it can be neglected. 
For lower q, however, it may become necessary to include a 
nontrivial transfer function Z. p if one aims for quantitative 
accuracy. It may be evaluated alternatively from eq. (|106[) . 



where a,/3 = 1,2 are spin indices and ei2 = —£21 = 1- 
Correspondingly, the propagator matrix reads 

detpVp 1 = (156) 

( h^>Qt a p , (-iq + |]g - cr)5 a p \ 

with 

-» 2 -^2 
det F = (iq + 7^- <r){-iqo + ^ - <r) + h%%. (157) 

The fermionic occupation number obtains as 



nF{q) = - : 



1 



tr 



, - 

8 a (3 , 



(158) 



For 4>q = we recover twice the result (|TT)) . corresponding 
to the two spin states. 

Let us evaluate n F {q) for T = 0. The zeros of detj? occur 
for 



XII. GAPS AND PSEUDOGAPS 

Besides a formal theoretical framework for systematic 
non-perturbative approximations, the exact flow equation 
for the occupation numbers can also be employed to mo- 
tivate comparable simply renormalization group improve- 
ments. Examples are the discussion in sect. IVII or the 
computation of the occupation numbers by a solution of 
the differential equation (|69p . with parameters evaluated 
at k 2 = q 2 . In this section we discuss gaps and pseudogaps 
in the light of such renormalization group improvements. 

Let us first discuss fcrmions in presence of an order pa- 
rameter or mean field. For definiteness, we consider a 
system of fermionic atoms, interacting with an effective 
bosonic di-atom fields <p. This field can represent micro- 
scopic molecules as well as Cooper pairs or other collective 
states with atom number two. For our purpose it is suffi- 
cient to consider a rather simple effective action for fermion 
fields tp with two "spin states" and <f> 



r 



A 

2M 



(153) 



The Yukawa or Feshbach coupling couples the fermions 
to bosons. For a constant real order parameter (j)(x) = 4>o 
this action is Gaussian and reads in momentum space 



1 



r = ^ / V (-q)V F (q)^(q) , i/>{q) 



f m 



(154) 



The inverse propagator Vf is a 4 x 4 matrix 



V F (q) 



(155) 



ql 



2M 



■Wo 



One finds 
n F (q) 



\2M 



\2M 



(159) 



-1/2 



(160) 

such that the Fermi surface is smoothened with a typical 
width set by the gap A = h^tpo. This generalizes to T 7^ 



with positive square root 



exp ( - ) + 1 




A 



(161) 



(162) 



In presence of a gap the fermion density n f does not vanish 
for T = and any finite a. For T = and a = the 
occupation numbers approach 



n F {q) = 1 



1 



\J q^ + AM 2 A 2 2MA 
According to cq. (|159[) the dispersion relation 
Hq)\ = e{q) 



(163) 



(164) 



is characterized by a minimal frequency set by the gap 

Mmin = A. 

Renormalization group improvement implies that for any 
q the gap A(q) should be evaluated at a characteristic mo- 
mentum scale k. For fermions this implies 



A(q) = h^Mk) , k 2 = \q 2 - 2Ma\. 



(165) 
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Effectively, the fermions with momenta distant from the 
Fermi surface, q = <jp + Sq , qp 2 = 2Ma, are influ- 
enced by an effective volume with linear size ~ where 
k 2 = \2qpSq + Sq 2 |. (The energy needed for a change 
Sq in momentum is k 2 /(2M), and k is a typical momen- 
tum of a boson whose absorption / emission can induce 
such a change.) The quantities h^ik) and 4>o(k) can be 
associated with the "running couplings" evaluated at the 
renormalization scale k. (This is precisely what happens if 
one investigates a fc-dependent particle density as in [§].) 
A non-vanishing gap A corresponds to 

A = h^(k = 0)4>a{k = 0) = h(k = 0)^/po (166) 

with renormalized Yukawa coupling h = h$l V~A. (We omit 
in our discussion the wave function renormalization for the 
fermions.) We observe that a gap does not necessarily re- 
quire a nonzero order parameter </> . A nonvanishing renor- 
malized order parameter ^fp~^ > is sufficient, provided the 
renormalized Yukawa coupling remains also nonvanishing 
for k 0. 

Furhermore, one may encounter situations where A(/c) 
differs from zero for some range k > ksn, while A(/c < 
ksn) = 0. For small enough ksn this is a typical situation 
for a pseudogap. Except for a very narrow region around 
the Fermi surface the occupation numbers np(q) behave as 
if a gap were present. Such a situation has been discussed 
in detail for antifcrromagnetic order in the two dimensional 
Hubbard model [l2j . A pseudogap occurs for temperatures 
between the critical temperature T c and the pseudocritical 
temperature T pc > T c . The presence of a pseudogap A(/c) 
simply reflects the existence of local order in domains with 
linear size ~ k^ 1 . 

Similar considerations hold for bosonic occupation num- 
bers. For non-relativistic bosons we may interprete A r = 
\/2Xpo as a gap for the radial mode, while the Goldstone 
boson remains gapless. More complex gaps may occur if 
different bosonic degrees of freedom play a role. A gap for 
the bosons also exists in the symmetric phase where the 
effective potential U(p) has a minimum at po = 0, with 
m 2 = dU/dp\ > 0. In this case one replaces in eq. (fT5)l 
— a — > m 2 . 



XIII. CONCLUSIONS 

In this paper we have concentrated on the more formal 
aspects of the functional integral approach to occupation 
numbers n(q) for interacting non-relativistic particles. Our 
study of bosons at zero temperature reveals pitfalls for 
computations beyond leading order perturbation theory - 
the Bogoliubov approximation. They are connected to the 
infrared and ultraviolet behavior of n(q). If insufficient care 
is taken, no meaningful result for the particle density n will 
be obtained. This is a very practical problem for precision 
computations for ultracold atoms, where the relation be- 
tween the particle density and the chemical potential a is a 
crucial quantity. One of the potential problems arises from 
an unphysical dependence of the relation between density 
and chemical potential, n{a), on the ultraviolet cutoff A. 



A second class of problems concerns the infrared behav- 
ior of the occupation numbers. For T = we have found 
that the behavior of n(q) for q — > differs qualitatively 
from the Bogoliubov theory for d = 1,2 and mildly for 
d = 3. This is due to the dominance of a term quadratic 
in the frequency in the full inverse propagator. Denoting 
with V the corresponding coupling we find 

n(q) = 4>l5(q) +n p (q), 

n p (q)~V~ 1/2 q- {1 - , i\ (167) 

The anomalous dimension r\ relates to the dense regime or 
Goldstone regime It is positive for d = 1 and vanishes 
for d = 2,3. As one of the consequences we find that 
the condensate depletion n p = J-n p (q) is not dominated 
by the extreme infrared modes with q — > 0. This allows 
for a clear conceptual separation between the condensate 
density n c = 0q and n p . For d — 3 the asymptotic behavior 
(|167[) sets in only for extremely small q, such that for most 
practical purposes the Bogoliubov approximation remains 
valid. 

For practical calculations with a good accuracy we 
propose to use the differential equation (JBSJ) for the q- 
dependence of n p {q). Correspondingly, n p can be com- 
puted from eq. ([70|) . The benefit of such a calculation 
is the absence of infrared or ultraviolet problems and the 
renormalization group improvement due to the use of run- 
ning couplings. We postpone a numerical solution of eqs. 
(gg) , (LTHJ) to future work. 

We have derived an exact functional renormalization 
group equation for scale-dependent occupation numbers. 
They can be employed to motivate a renormalization group 
improvement in simpler settings. The effective couplings 
appearing in the evaluation of the occupation number n(q) 
for a particular momentum q correspond to running cou- 
plings, evaluated at a renormalization scale k = q (with 
modifications in presence of a Fermi surface). This allows 
for an easy description of pseudogaps. 

APPENDIX A: One loop calculation of Z$ 

We define the wave function renormalization Z$(u},q) by 
the w-derivative of the 2 — 1-element of the inverse propa- 
gator 

Z 4> {u,q) = -i^- j P 21 . (A.l) 

For nonzero <po the one loop contribution to P21 has a con- 
tribution from effective cubic couplings ~ Xcf>o 

AP 21 (u;,q) = -\ 2 ^ 2 f f 
J Pa Jp 

|3Gii(po,p)(G2i(po + iu,p + q) + Gi 2 (po - iw,p - qf) 
+G 22 (po,p)(Gi2(po + iu,p + q) + G 2 i(po - iu,p- q))\ 

(det- 1 {p+q)+det- 1 { P -q)), (A.2) 
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with 



det(p + q) 



AM 2 



(P + q) 2 [(P + q? + 4A/A 2 ] + ( po + iu) f, 

(A.3) 



We first consider the special case q = 



i{»[j(p» + J(f,-»)]}, (A.4) 



where 



J(P,oj) = / det ^det \p + q) 

Jpn 

(j?l+B)- x ((p +iuf+B)-\ 



(A.5) 



The integrand of the po-integration has simple poles at 
Po = ±i\/~B and po = -k.i\f~B — iu> and we take uj > 
without loss of generality. The first pole in the upper half 
plane at po.i = iVB has residuum 



2^B~uj(2s[B + uj) ' 



(A.6) 



whereas a second pole in the upper half plane at po,2 
i(>/~B~ — uj) exists only for oj < y/B, with residuum at 



1*2 



2VBuj(2VB - uj) 



(A.7) 



This yields 



J(lo > 0) = — 



2-^Bui V 2-/B + uj 2\[~B - uj 



or, for arbitrary uj, 

K{uj) = uj(J(u) + J{-u)) 



{VB-uj) 
(A. 



VB(4,B-uj 2 ) 

sign(ij) 
s/B{2VB+\uj\) 



for U < V^B 



for U > VB 



(A.9) 



Both if and dK/duj are discontinuous at \u)\ = \/B. 

At first sight the situation seems to simplify for uj = 
where 



1 



dK 



dw \uj=o 2B\fB^ 



(A.10) 



resulting in 



With cr = 4MA^0q we have to solve the remaining integral 
over y = \p\ 



3,(0,0) = 4^Af 2 A^ 2 /(a), 

oo 

= / ^/- 4 Uy 2 + or 112 - y (y 2 + °r m 



3 , \ d—\ _ 

-d-5)a— F d . 



DC 



(A.12) 



This yields 

Z+(p, 0) = -2 d - 3 (10 - id)v d F d (\^M)i 4~ 2 . (A.13) 

However, the integral Fd is infrared finite only for d > 3 
and ultraviolet finite for d < 4. For the realistic cases with 
d < 3 the limit w — > has to be done more carefully. Even 
for very small \uj\ eq. (| A. 10|) can only be used for \[B > \u>\ 



z > 



2M\u)\ 



a 



(A.14) 



This modifies effectively the factor Fd in eq. (|A.13[) , which 
now becomes a function of uj/a. The contribution to Fd 
from this range is approximately 



F + - 



1 

3-d 



2Mu 



,d-3 



hi 



2JU* 



for d < 3 
for d = 3 



(A.15) 



and diverges for uj — ► 0. In consequence, one has Z^{u>, 0) ~ 

3 _ 

A^0o for d < 3 with logarithmic corrections for d = 3. The 
range z < 2Mui/a (for > 0) gives an infrared finite 
contribution for d > 1 since for z — > one has dK/duj « 
l/(w 2 \fB) and the r.h.s. of eq. (jA.llj) is multiplied by a 
factor 2B/w 2 z 2 cr 2 /(2Af 2 L i ; 2 ). For d = 1 the one loop 
expression for Z^(u), 0) remains infrared divergent even for 

> 0. In any case, the negative one loop correction (|A.13|) 
drives Z^uj, 0) rapidly towards zero as uj decreases, making 
the use of renormalization group improvement mandatory. 
This has been done in Q. 

Let us now turn to the evaluation of Z^iio^q) for general 
q and uj. By suitable momentum shifts eq. (|A.2[) can be 
written as 



AP 21 



iuj\^ 



H = 



M 



det 



2M + £- 



H. 



det" 



(A.16) 
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where 
det 



B+ = 



B± 
1 



Po ± ~ w 



(A.17) 



4M 



2 



For the ^-integration we need the poles which are situated 
at p ,i = +iyfB^ - \u , p oa = +i^/~B~ + |w, p ,3 = 
—iy/B+ — |w , po,4 = —iy/B- + ^lj. For |w| < 2B ± one 
finds 



H = i(n + r 2 ), 

in = {2 v ^b7(b_ -(V^ 

ir 2 = {2v^Br(s + - (v/fiT + w) 2 )} . (A.18) 



For g = one recovers the formulae discussed above. 
For g = |g| — > oo we can neglect u such that 



B± 



1 



4M 2 



(A.19) 



We may evaluate eq. (|A. 18|) for uj = 
r-i 



2V-B+B_( J B + +B_) 

2X 2 



1 



2V2A/ 3 



P 2 + ^) -l^/r 



(A.20) 



This yields 



AZ& = 2V2A 2 2 M 2 



^o 1 

2N 2 



m 2 



-1/2 



(A.21) 



The p-integral is finite for d < 4 and dominated by values 
p 2 pa <f 2 /4. As a result it is proportional Udg d ~ 4 /(<i— 4) and 
we find an expression as in eq. (|4"T|) . Solving the integral 
explicitely determines the coefficient 7 in this equation. 



We denote the inverse boson propagator by P^ — G^ 1 . 
We consider here an approximation where is given by 
a classical part and a part AP^ induced by the fermion 
fluctuations 



P^ = iq + ^-2a + AP^. 



(B.l) 



(The boson mass is twice the fermion mass M and the 
chemical potential twice the fermion chemical potential. 
This corresponds to di-atom molecules.) This inverse prop- 
agator has been extensively studied in [§]. We are inter- 
ested here in its value for analytically continued momenta 
go = «w and at T — 0. We want to show here that the 
fermion fluctuation contribution to the boson (or molecule) 
occupation number riM{q) vanishes in absence of a gap 
A = h$(f>o = h</,(t>o (for T = 0). 
Using the shorthands 



uj = M(ui + 2a), q = 2Mq' , q = Mq 



(B.2) 



one finds 
AP^iqa + iuj,q) 



hlM 



<7 ' 



dq 



-\ 2 

q - 7j ) + l io - u + m 



q I 

9 + o ' +m-^-iq Q 



(uj^0,q^0,q o ^0) 



(B.3) 



The go integration has poles at 
2 



9o = « 



qo , -i 



(B.4) 

and can be performed by closing the integral either in the 
upper or lower half plane 



AP d> = -hiM 



g ' 



q u + — -u + tq 



[t-2Ma]~ 1 \. 



(B.5) 



We concentrate here on d — 3 and a < 



APPENDIX B: Fermion boson coupling 

In this appendix we briefly discuss the effects of fermion 
fluctuations on the bosonic occupation numbers in a cou- 
pled system as given by eq. (|153[) . For this purpose the 
effective action has to be extended by a piece quadratic 
in the boson field <f> and by suitable boson interactions. 



AP rh = 



A 2 



K M 

4tt 2 ./ \(q' - iA)(q' + iA) 
2Ma 



(q> - iy/-2Ma)(q / + W-2Ma) 

with A 2 = q 2 /4 — uj + iq . We choose A to be the root 
with positive real part. (This always exists for go 7^ 0.) 
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The q'-integration yields 

hlM 3 / 2 



4tt 



/ ->2 \ 1/2 



(B.6) 



1/2 

where ( ) / indicates the root with positive real part. We 
therefore obtain the propagator in the molecule phase (a < 
0) 



4M 



(B.7) 



1/2 



We want to show that the bosonic occupation number 
vanishes 



™a/(<z) = / G(q ,q) = 0. 



(B- 



For this purpose we consider dnM{q) / 'dti 2 ^. For this quan- 
tity we can close the integral on the lower half plane. 



The integral is analytic in the lower half plane, implying 
driM(q 2 )/ 'dh^ — 0. Since for = the theory is free we 
infer for all 2 > that riM(q; = 0) = 0. Therefore 
ntaiq) vanishes for all h 2 ^. 

For the atom phase one has a = and the bosonic prop- 
agator becomes 



G(q ,q) = [iqo + jt? +rrr 



h 2 M 3 / 2 

47T 



2 \ I/ 2 ' 



AM J 



(B.9) 



The argument showing Um® = proceeds similar as 
above. 

In the limit (jB.ip the boson occupation number vanishes 
for T = as long as there is no gap. A non- vanishing 
fermion contribution to the boson numbers only arises from 
modifications of the fermion propagator in eq. (|B.3j) due 
to a non- vanishing gap A = h^o- Of course, the fermion 
fluctuation contribution to riM{q) is nonzero for T > 
even in absence of a gap |8|. Furthermore, one has the 
contribution from boson fluctuations as discussed in the 
main part of this paper. 
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